Spatial correlations of mobility and immobility in a glassforming Lennard-Jones liquid 
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Using extensive molecular dynamics simulations of an equilibrium, glass-forming Lennard-Jones 
mixture, we characterize in detail the local atomic motions. We show that spatial correlations exist 
among particles undergoing extremely large ("mobile") or extremely small ("immobile") displace- 
ments over a suitably chosen time interval. The immobile particles form the cores of relatively 
compact clusters, while the mobile particles move cooperatively and form quasi-one-dimensional, 
string-like clusters. The strength and length scale of the correlations between mobile particles are 
found to grow strongly with decreasing temperature, and the mean cluster size appears to diverge 
at the mode-coupling critical temperature. We show that these correlations in the particle displace- 
ments are related to equilibrium fluctuations in the local potential energy and local composition. 

PACS numbers: 02.70.Ns, 61.20.Lc, 61.43.Fs 



I. INTRODUCTION 

The bulk dynamical properties of many cold, dense 
liquids differ dramatically from what might be expected 
from extrapolation of their high temperature behavior 
0. For example, many liquids cooled below their melt- 
ing temperature exhibit rapid non-Arrhenius increases of 
viscosity and relaxation time with decreasing tempera- 
ture, and two-step, stretched exponential decay of the 
intermediate scattering function F(q, t). Such behavior 
is often discussed as a "signature" of the approach to the 
glass transition. It has long been a central goal of theories 
of the glass transition to account for these bulk phenom- 
ena in terms of the microscopic dynamical motions of 
the molecules of the liquid. As a consequence, computer 
simulations of supercooled liquids, in which this micro- 
scopic information is immediately available, are increas- 
ingly used to complement theoretical and experimental 
efforts. In particular, simulations in recent years have 
focused on the study of "dynamical heterogeneity" to 
understand the microscopic origin of slow dynamics and 
stretched exponential relaxation in glass-forming liquids 



Recently we reported the observation of dynamical het- 
erogeneity and also cooperative molecular motion ||] 
in extensive molecular dynamics simulations of a super- 
cooled Lennard-Jones (LJ) mixture. These spatially cor- 
related dynamics were observed in a regime of tempera- 
ture T, density p and pressure P for T above the dynam- 
ical critical temperature T c obtained PJlCj ] from fits by 
the ideal mode coupling theory (MCT) Jll[| to other data 
on the same system. The principle goals of the present 
paper are twofold: (1) To test directly for spatial corre- 
lations of particles assigned (according to their displace- 



ment over some time) to subsets of extreme mobility or 
immobility, and (2) to establish connections between this 
"dynamical heterogeneity" and local structure. 

This paper is organized as follows. In Section || we 
present relevant background information, and in Sec- 
tion III we describe the model and computer simulation 



techniques. 



In Sec. IV, we examine the bulk dynamics 

In Sec. 



and equilibrium structure of the liquid, in Sec. |Vj we 
examine the mean square displacement and analyze the 
shape of the time-dependent distribution of particle dis- 
placements to define a time scale which we use to study 
dynamical heterogeneity throughout the remainder of the 



paper. In Sec. VI we group particles into subsets accord- 



ing to the maximum displacement they achieve on the 
time scale defined in the previous section, and show that 
particles of extremely high or low displacement are spa- 
tially correlated. From this we are able to identi fy a 
length scale that grows with decreasing T. In Sec. VII 



we show that fluctuations of the local mobility are cor- 
related to fluctuations of the potential energy, or alter- 
natively to fluctuations in the local composition of the 
liquid. In Sec. VIII we examine certain time dependent 
quantities associated with the observed dynamical het- 
erogeneity, and finally in Sec. IX we conclude with a dis- 
cussion. 



II. BACKGROUND 

It has been proposed that the stretched exponential 
behavior exhibited by the long time relaxation of F(q, t) 
can be attributed to a sum of many independent local ex- 
ponential relaxations with different time constants, i.e., 
to a distribution of relaxation times . This interpre- 
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tation is one form of the so-called "heterogenous" sce- 
nario for relaxation [|6|,^2|-p7t . A number of recent experi- 
ments [p~3| have shown that in liquids such as orthoter- 
phenyl and polystyrene within 10 K of their glass tran- 
sition temperature T g , subsets of molecules rotate slowly 
relative to the rest of the molecules on time scales long 
compared with collision times, but shorter than the re- 
laxation time of density fluctuations. These liquids were 
thus termed "dynamically heterogenous." None of these 
experiments were able to explicitly demonstrate whether 
slow molecules are spatially correlated, but typical dis- 
tances over which slow molecules may be correlated were 
inferred [fl3| . 

There have been numerous attempts to indirectly mea- 
sure a characteristic length scale over which molecular 
motions are correlated at the glass transition both in ex- 
periments @-H and in simulations Donth jlf] 
relates the distribution of relaxation times in systems ap- 
proaching their glass transition to equilibrium thermo- 
dynamic fluctuations having a characteristic size of ~ 3 
nm at T a . Thermodynamic measurements on orthoter- 
phenyl , and dielectric measurements on salol |^0| , N- 
methyl-e-caprolactan and propylene glycol pl| , showed 
a shift in T g due to confinement in pores of the order 
of a few nanometers. Mountain showed that the size 
of regions that support shear stress in a simulation of a 
glass-forming mixture of soft spheres grows with decreas- 
ing temperatures. Monte Carlo simulations of polymer 
chains in two dimensions demonstrated strong finite size 
effects on diffusion |22j. A number of experiments and 
simulations on polymers confined to thin films all found 
a shift of T g due to confinement ]23|-|27|] . These effects 
have all been attributed to the presence of cooperatively 
rearranging regions that grow with decreasing T. How- 
ever, the origin of this characteristic length has never 
been shown explicitly. In particular, the connection of 
the characteristic length to a cooperative mechanism of 
molecular motion has not been experimentally demon- 
strated. 

The intuitively-appealing picture of cooperative molec- 
ular motion was proposed in 1965 by Adam and 
Gibbs |p8|| . In their classic paper, they proposed that 
significant molecular motion in a cold, dense fluid can 
only occur if the molecules rearrange their positions in a 
concerted, cooperative manner. They postulated that a 
glass-forming liquid can be viewed as a collection of in- 
dependently relaxing subvolumes within which the mo- 
tion of the particles is cooperative. As the temperature 
of the liquid is lowered, the number of particles involved 
in cooperative rearrangements increases. If structural re- 
laxation occurs through the cooperative rearrangement of 
groups of molecules, the liquid observed over a time-scale 
shorter than the structural relaxation time will appear as 
a collection of regions of varying mobility. These predic- 
tions can be tested by selecting subsets of molecules that 
relax slower (or faster) than the average, and determining 



whether the molecules in a subset are randomly scattered 
through the sample or tend to cluster in a characteristic 
way. 

The explicit connection between dynamical hetero- 
geneity and cooperative motion is only recently being in- 
vestigated experimentally in detail p0[ ]. However, there 
have been a number of recent computational investiga- 
tions addressing these issues. For example, Muranaka 
and Hiwatari showed that displacements of particles 
measured over a timescale of the order of 5 collision 
times are correlated within a range of about two inter- 
particle distances in a two-dimensional binary mixture 
of soft disks below the freezing point. Wahnstrom 0] 
showed that hopping processes in a strongly supercooled 
binary mixture are cooperative in nature. Hurley and 
Harrowell B identified fluctuating local mobilities in a 
supercooled two-dimensional (2-d) soft-disk system, and 
showed an example of correlated particle motion on a 
timescale of the order of 20 collision times. Mountain 
demonstrated similar correlated particle motion in a 
2-d supercooled Lennard-Jones mixture. By examining 
the time at which two neighboring particles move apart 
in 2-d and 3-d simulations of a supercooled soft-sphere 
mixture, Yamamoto and Onuki demonstrated the growth 
of correlated regions of activity || . They further stud- 
ied the effect of shear on these regions Q, and showed 
that the size of the regions diminished in high shear. The 
clusters of "broken bonds" (denoting pairs of neighboring 
particles that separate beyond the nearest neighbor dis- 
tance) identified in that work are similar in some respects 
to the clusters of highly mobile particles in a 3-d binary 
Lennard-Jones liquid reported previously by us [Q, and 
described in detail in the present paper. The connec- 
tion between the clusters of Ref . J7| , which demonstrate 
a form of dynamical heterogeneity, and cooperative par- 
ticle motion, was shown in Ref. M. 

III. SIMULATION DETAILS 

We performed equilibrium molecular dynamics (MD) 
simulations of a binary mixture (80:20) of N = 8000 
particles in three dimensions. The simulations were per- 
formed using the LAMMPS molecular dynamics code ]3C| l 
which was designed for use on distributed memory par- 
allel machines. LAMMPS partitions particles (atoms or 
molecules) across processors via a spatial decomposition 
|nj whereby each processor temporarily "owns" parti- 
cles in a small fixed region of the simulation box. Each 
processor computes the motion of its particles and ex- 
changes information with neighboring processors to com- 
pute forces and allow particles to migrate to new proces- 
sors as needed. 

The 6400 particles of type A and 1600 particles of type 
B interact via a 6-12 Lennard-Jones potential, 
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V a p(r) = 4e Q/3 



in 



where aj3 € {A, B}. The interaction forces between par- 
ticles are zero for all r > r c = 2. ,5a a a- Both types 
of particles are taken to have the same mass to. The 
Lennard-Jones interaction parameters e aj( g and a a a for 
this mixture are: €aa = 1.0, cab = 1.5, £_bb = 0.5, 
oaa = 1-0, (Jab = 0.8, (t bb = 0.88. Both the relative 
concentration of particle types and the interaction pa- 
rameters were chosen to prevent demixing and crystal- 
lization ||. Throughout this paper, lengths are defined 
in units of ctaa, temperature T in units of eAA/ks, and 
time t in units of \J a AA m/tAA- 

The simulations for each state point (P, T, p) are per- 
formed in three stages. First, a constant NPT adjust- 
ment run is performed by coupling the system to stochas- 
tic heat and pressure baths to bring the system from a 
nearby state point (usually the previously simulated state 
point) to the desired state point ]32] ]. Second, a constant 
NVT equilibration run is performed to test for unwanted 
drifts in pressure P or potential energy U [ |33| . If no drift 
is observed, the final state of the system is considered to 
represent an equilibrium state of the system. Third, a 
constant NVE data-gathering run is performed using the 
final equilibrated state obtained from the second stage, 
and snapshots containing the particle coordinates and 
velocities are taken at logarithmic time intervals during 
the run. In this stage, the equations of motion are in- 
tegrated using the velocity Verlet algorithm with a step 
size of 0.0015 at the highest temperature, and 0.003 at 
all other temperatures. All quantities presented here are 
calculated in this third stage. The analysis is performed 
by post-processing the snapshot files, which number as 
many as several thousand for the lower temperatures. 

We simulated nine state points along a path in P, T, p 
that is linear when projected in the (P, T) plane. This 
path was chosen so that we would approach, from high 
temperature, the mode-coupling critical point T c — 
0.435, P c = 3.03, p c — 1.2 || along a path different from 
that used in Ref . || . Table 1 shows the values of P, T, 
and p for each state point. 

For state points far above T c (e.g. runs 1-5), the data- 
gathering runs required more cpu time than the equilibra- 
tion. For state points nearer T c , the equilibration stage 
was the most time consuming. In these cases the NVT 
stage of the simulations showed a slight drift of the pres- 
sure over time. To shorten the time required for complete 
equilibration, we estimated a new volume or temperature 
to create a nearby state point that we expected would be 
very nearly equilibrated with the current particle config- 
uration. Then we instantaneously scaled the positions or 
velocities of the system (ie. adjusted the volume or tem- 
perature) and began another constant NVT run to test 
for equilibration. By iterating this procedure a few times 
we were able to find an equilibrated state point within 



0.03% of the desired P and T at low temperatures. 

At the lowest T studied (T = 0.4510), the total run 
time following equilibration is 1.2 x 10 4 time units. Thus, 
assuming Argon values for the parameters in Eq. 3.1, the 
data presented here extend up to 25.8 ns. 



IV. STRUCTURE AND BULK RELAXATION 

In this section, we show that the simulated liquid 
exhibits the characteristic features of an atomic glass- 
forming liquid. 

Structural relaxation may be probed experimentally 
by the intermediate scattering function P(q, t), which is 
both the spatial Fourier transform of the van Hove corre- 
lation function G(r, t) and the inverse time transform of 
the dynamic structure factor S*(q, ui) |34f| . In a computer 
simulation, the self (incoherent) part of the intermediate 
scattering function P s (q, t) may be calculated directly 
from 



P s (q,*) 



(2) 



where rj(t) is the position of particle j at time t, and 
(■ • •} indicates an average over independent configura- 
tions. This quantity describes the relaxation of density 
fluctuations due to single particle displacements on an 
inverse length scale 2n/q, where q = |q|. If we assume 
rotational invariance of the system, F s (q, t) depends only 
on q. The time dependence of F s (q, t) for the A particles 
for q = q max is shown in Fig. [lj (Throughout this paper, 
q is chosen as q m ax, the position of the first maximum of 
the static structure factor S(q,0)). At high T, F s (q,t) 
decays to zero exponentially. As the system is cooled, 
F s (q,t) develops a plateau that separates a short time 
relaxation process from a long time relaxation process. 
This plateau indicates a transient "localization" of par- 
ticles in the "cages" formed by their neighbors, and is a 
characteristic feature of all glassforming liquids. 

The mode-coupling theory developed for supercooled 
liquids by Gotze and Sjogren makes a number of predic- 
tions concerning the decay of the intermediate scattering 
function O). These predictions have been tested and 
verified for the LJ potential used here in a regime of P, 
T, and p similar but not identical to that simulated here 
There it was shown, e.g., that the early and late 
/3-relaxation regimes are well described by power laws, 
and that the late time behavior of F s (q, t) exhibits time- 
temperature superposition with a time constant T a that 
diverges as a power law as T approaches T c ~ 0.432, with 
exponent 7 ~ 2.7. The diffusion constant was found to 
scale as D ~ (T — T c )~ 7 , with 7 = 2.0 for the A particles, 
7 = 1.7 for the B particles, and T c = 0.435. 

The simulations performed in the present work extend 
from a point in the phase diagram where two-step re- 
laxation begins to emerge, down to a state point that 
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is within approximately 4% of T c . Over this range, we 
find that r Q increases by 2.4 orders of magnitude, and 
fits well to the power law form found in Ref. 0, with 
approximately the same critical temperature and critical 
exponent. 

It is well known that although relaxation becomes 
strongly nonexponential and relaxation times increase by 
many orders of magnitude as a supercooled liquid ap- 
proaches a glass transition, changes in the static structure 
of most liquids are far less remarkable. To demonstrate 
this for our system, we examine the pair correlation func- 
tions g a p{r) given by 



V 



for a (3 and 



N a Np 



(3) 



je/3 



1,3 

where N a (Np) is the total number of particles of species 
a (f3). With this normalization, g a p(T) converges to unity 
for v — > oo in the absence of long range correlations. 
Assuming rotational invariance, the correlation functions 
do not depend on the direction of the vector r, but only 
on the distance r = |r|. 

In Figs. |,| and g we show the pair correlation func- 
tions gAA{f), gAB(r), and gss(r) for three temperatures. 
The figures show that these functions do not change dra- 
matically as a function of the state point. As the tem- 
perature is lowered, the main effect on all three func- 
tions is that the maxima and the minima become slightly 
more pronounced. Additionally, the second maximum of 
gAA{r) and sus(f) at low T shows a splitting that has 
commonly been interpreted as a signature of an amor- 
phous solid, although at these state points our system 
is an equilibrium liquid. Recently, evidence has been re- 
ported [[35| that in a 2-d system of hard-disks the splitting 
of the second peak in the pair correlation function is due 
to the formation of regions with hexagonal close-packed 
order. 



V. SINGLE PARTICLE DYNAMICS 

Having established that the model liquid studied here 
exhibits the characteristic bulk phenomena of a glass- 
forming liquid, we examine in this section the distribu- 
tion of individual particle motions. 

The most basic dynamical bulk quantity that is easily 
accessible to simulation is the particle mean square dis- 
placement (MSD), (r 2 (t)). Because we are investigating 
a binary mixture, we refer in the following to a MSD for 
the A particles and a MSD for the B particles. At high T, 
the MSD for both species exhibits two distinct regimes 



(see Fig. ||). In the short time limit (regime I) the MSD 
is ballistic, i.e. (r 2 (t)) oc t 2 . For longer times (regime 
III), the MSD is diffusive, i.e. (r 2 (t)) oc t. As the system 
is cooled, an intermediate regime (II) between these two 
limiting behaviors develops. Before entering the diffu- 
sive regime, (r 2 (t)) exhibits a plateau, analogous to the 
plateau in the intermediate scattering function, that like- 
wise arises from a transient "caging" of each particle by 
its neighbors. As seen in the figure, the time the sys- 
tem spends in the plateau depends strongly on T, and 
increases with decreasing T. The MSD for the B par- 
ticles (not shown) exhibits qualitatively the same time 
dependence as shown in Fig. 5 , but the diffusive regime 
is reached at shorter times, and the diffusion constant is 
larger, than for the A particles [0. This difference can 
be explained by the different sizes of the A and B parti- 
cles and by the fact that the interaction constant cbb is 
smaller than caa- 

In this paper, we are interested in whether spatial cor- 
relations exist between particles that exhibit either ex- 
tremely large or extremely small displacements over some 
time interval. To determine this, we must first define the 
time interval over which the particle displacements will 
be monitored. Obviously, displacements may be mon- 
itored over any time interval, from the ballistic regime 
to the diffusive regime. To see whether there is a natu- 
ral time scale on which the particle displacements might 
exhibit a particularly strong correlation, we turn to the 
self part of the van Hove correlation function, G s (r,t), 
which gives the probability to find a particle at time t at 
a distance r from its position at t — J34|: 



G s (r,t) 



^(r + r 4 (0)-r,(i)) 



(5) 



Due to the rotational symmetry of the system, G s (r,t) 
is a function of the modulus r of the vector displacement 
r. The quantity Amr 2 G s (r, t) , which gives the number of 
particles located a distance r from their original position 
at time i, is shown in Fig. [(] for the A particles for three 
different times at the lowest T. Also shown in Fi 
the Gaussian approximation 4 7rr 2 G°(r, i) , where 




G°(r,t) = 



2ir{r 2 (t)) 



exp 



3r 2 



2(r 2 (t)) 



(0) 



and where (r 2 (t)) is equal to the measured one. The 
Gaussian form appears to be a good approximation to 
G s (r,t) at both short and long times. However, it is 
apparent from the figure that G s (r, t) is significantly dif- 
ferent from G°(r,t) at intermediate times. In particu- 
lar, while many of the particles have traveled less than 
would be expected from the knowledge of (r 2 (t)) alone, a 
small number of particles have traveled significantly far- 
ther. As a result, at intermediate times G s (r,t) displays 
a long tail that extends beyond one interparticle distance 
at T= 0.4510 (cf. Fig. g). 
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This "long tail" behavior is most pronounced at a 
time f * when G s (r, f ) deviates most from a Gaussian (cf. 
Fig. Jr| as characterized by the "non-Gaussian" parame- 
ter pQl, 



a 2 (t) 



3(r 4 ft)) 
5<r 2 (f)} 2 



-1, 



3. 



(7) 



While one can define many different non-Gaussian pa- 
rameters, this particular one involves the lowest possible 
moments. With this definition, 02(f) is zero if G s (r,t) is 
Gaussian. If a distribution has a tail that extends to dis- 
tances exceeding those for a Gaussian distribution with 
the same second moment, all higher order moments of 
this distribution will exceed those of the corresponding- 
Gaussian, and consequently 02(f) will assume positive 
values. In Fig. || we show 02(f) for various T for the A 
particles. As expected, 02(f) is zero at short times, then 
becomes positive, exhibits a maximum, and finally goes 
to zero at long times. As T decreases, the position of the 
maximum f* shifts towards longer times, and the height 
of the maximum increases. For all T, we find that 
f* corresponds to times in the late-/3/early-o relaxation 
regime. Furthermore, by dividing 02(f) by o|, and divid- 
ing f by f * , one can show |37j that all curves collapse onto 
a single master curve for all times larger than the micro- 
scopic time (where they already collapse before scaling) 
P, ^0|j3^ ]. This data collapse, which is likely related to 
the time-temperature superposition exhibited by the in- 
termediate scattering function (but not trivially related, 
since f* does not scale linearly with t q ), suggests that f* 
is in some sense a characteristic time for this system [[37) . 
Note that f * is orders of magnitude larger than the micro- 
scopic "collision time" |39| r; for example, at T = 0.4510 
f* = 155.5 and r = 0.09. 

We see from an analysis of the temperature-dependent 
distribution of particle displacements at various times, 
and the calculation of the time-dependence of the non- 
Gaussian parameter, that the single particle dynamics is 
most non-Gaussian — and displays the widest range of 
possible behaviors — on the timescale f*. The interval 
from zero to f* thus provides a convenient choice over 
which to monitor the particle displacements and study 
their correlations because (i) since f* is the time at which 
the distribution of particle displacements is broadest, it 
may also be when the liquid is likely to be most "dy- 
namically heterogeneous" ; and (ii) f * is well-defined and 
easily calculated. Thus, throughout this paper we will 
use the time window from zero to f* as the time interval 
over which the particle displacements are calculated, and 
over which we investigate dynamical heterogeneity. 



VI. ANALYSIS OF SPATIAL CORRELATIONS 
OF PARTICLE DISPLACEMENTS 

In this paper, we are interested in studying the extreme 
behavior of the individual particle motion, from the ex- 
tremely mobile to the extremely immobile. In Refs. |?j 
and ||, we defined the magnitude of the displacement 
itj(f, f*) = |r.j (t* + f ) — r,(i)| of particle i in a time inter- 
val f * , starting from its position at an arbitrarily chosen 
time origin f , as a measure of the mobility of the i-th par- 
ticle. At t* , the distribution of Ui values is given by the 
self part of the van Hove correlation function, G s (r, f*), 
where r = u (cf. Fig. 0). In Refs. Q|| a subset of 
"mobile" particles was defined by selecting all the par- 
ticles that in the interval [0,f*] had traveled beyond the 
distance r* where G s (r, f*) exceeds G°(r, f*). With this 
definition, "mobile" particles are those that contribute 
to the long tail of the van Hove distribution function at 
the time f* (cf. Fig. @). In Refs. @,||, it was shown that 
mobile particles selected according to this rule tend to 
cluster [Q, and move cooperatively This definition 
of mobility given by the magnitude of particle displace- 
ment is thus sufficient to establish the phenomena of both 
dynamical heterogeneity and cooperative motion. 

Intuitively, we think of immobile particles as those par- 
ticles which are trapped in cages formed by their neigh- 
bors. Nevertheless, particles do not sit at one position; 
they essentially "oscillate" back and forth within the cage 
formed by their neighbors. To study correlations between 
the most immobile particles, we need a definition of mo- 
bility which allows us to select the particles for which the 
amplitude of this oscillation (ie. the maximum displace- 
ment of the particle) is the smallest. In this paper, we 
therefore define the mobility /i^ (f ) of the i-th A particle 
as the maximum distance reached by that particle in the 
time interval [f, t + t*]: 



Mi(f) 



V e [o, t .]{K-(f' + i)-ri(f)|} 



(8) 



This new definition of mobility, which we use throughout 
this paper, allows us to examine different subsets of par- 
ticles, from the most to the least mobile, in the same way. 
As a compromise between examining the most extreme 
behavior and including enough particles to obtain good 
statistics when examining their spatial correlation (i.e. 
maximizing the signal-to-noise ratio), we will examine 
the 5% most mobile and 5% least mobile particles. Thus 
we define as "mobile" the 5% of all particles having the 
highest values of //(f), and "immobile" the 5% having the 
lowest value. Note that this new definition of mobility 
does not qualitatively change the results obtained pre- 
viously in |@,||, provided that the new definition selects 
approximately the same fraction of the sample as the def- 
inition previously used (approximately 5.5% in Ref. 0) 
Compare, for instance, Figs, [ll] and 12 with Fig. 3 of 
Ref. §. 
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The subsets of mobile particles selected using the defi- 
nition of Ref. and that used here have a large overlap, 
since particles that have moved relatively far at some 
time in the interval [0, £*] are likely to remain relatively 
far at the end of the interval. However, subsets of im- 
mobile particles selected with the two different rules do 
not have as large an overlap, since a particle with a small 
displacement at some time may have previously traveled 
far, and then returned to its original position. The dis- 
tribution 4irfi 2 P(fj,, t*) at t* is shown in Fig. M. For com- 
parison, the probability distribution Airr 2 G s (r, t*) is also 
shown. Note that, although at t* particles can be found 
arbitrarily close to their position at t = 0, P{[i : £*) is zero 
for n < 0.17. 

In Fig. we show the 320 mobile particles (light 
spheres) and the 320 immobile particles (dark spheres) 
at the beginning of an arbitrary time interval [t, t + t*] 
for one configuration at T = 0.4510. The other 7360 
particles are not shown. The figure shows that particles 
of similar mobility are spatially correlated and that par- 
ticles with different mobility tend to be anticorrelated. 
These correlations can be quantitatively studied by cal- 
culating static pair correlation functions between parti- 
cles belonging to the different subsets. 

In Fig. [ll] we show the pair correlation function 
<?mm(?") between mobile particles for four different tem- 
peratures. gMAi(r) is defined by Eq. [I] with the sum 
restricted to the mobile particles.. For all T, <?mm(?") is 
appreciably higher that the average gAA{r) (cf. Fig. ||) 
for all r. The "excess" correlation given by the ratio 
T( r ) = [gM m (r) / gAA (?*)] — 1 is plotted as a function of 
r in Fig. |l2|. With the exception of the excluded volume 
sphere of the LJ potential, T(r) > at intermediate dis- 
tances and converges to zero for large r. It is clear from 
the figure that the total excess correlation, given by the 
area under the curve, increases with decreasing T. 

We can obtain an estimate of the typical distance over 
which mobile particles are correlated by identifying clus- 
ters of nearest-neighbor mobile particles jh|. To do this, 
we use the following rule: two particles belong to the 
same cluster if their distance at t = is less than r nn , 
the radius of the nearest neighbor shell, which is defined 
by the first minimum in gAA{f) and has a weak temper- 
ature dependence. In our hottest run r nn = 1.45, while 
in the coldest run r nn — 1.40. The distribution P(n) of 
clusters of size n is shown in Fig. Although most of 
the clusters have only a modest size, the data show that a 
significant fraction of the mobile particles, which them- 
selves make up only 5% of the sample (320 particles), 
are part of big clusters. For instance, at T = 0.4510, 
there is typically at least one cluster in each configuration 
that contains s=s 100 particles. For that T, P(n) ~ rT r 
with t — 1.86. In the inset we show the mean cluster 
size S — n2 P{ n )/Yl nP(n) Hj, plotted log-log versus 
T — T c , where T c = 0.435 is the fitted critical temperature 
of the mode coupling theory djnj . Although there is less 



than a decade on either axis, the figure shows that the 
temperature dependence of S is consistent with a diver- 
gence at T c of the form S ~ (T - T c )-f, with 7 « 0.618. 
Note that MCT makes no predictions about clustering or 
the divergence of any length scales as the critical point 
is approached p2[ . 

To test the sensitivity of the apparent percolation tran- 
sition at the mode-coupling temperature, we repeat the 
cluster size distribution analysis for the 3% and 7% most 
mobile particles. For each subset, the mean cluster size 
S is shown vs. T - T c in Fig. [[§ The best fit of 
S ~ (T - T p )-~< to each set of data gives T p = 0.440 
for the set containing the 3% most mobile particles, 
T p = 0.431 for the set containing the 5% most mobile 
particles, and T p = 0.428 for the set containing the 7% 
most mobile particles. However, within the accuracy of 
the data the three sets are also consistent with a diver- 
gence at T c . If we further increase the fraction of mobile 
particles beyond the fraction corresponding to a random 
close- packed percolation transition [ 13| , the mobile parti- 
cles percolate and most of the mobile particles are found 
in a single cluster that spans the whole simulation box. 

In Fig. |l5| we show one of the largest clusters of mobile 
particles found in our coldest simulation. It is evident 
from the figure that these clusters cannot be described 
as compact, as often supposed either implicitly or explic- 
itly in phenomcnological models of dynamically hetero- 
geneous liquids Instead, the clusters formed by 
the mobile particles appear to have a disperse, string-like 
nature. As discussed in Q, a preliminary calculation of 
the fractal dimension of the clusters, although hampered 
by a lack of statistics, indicates that the clusters have a 
fractal dimension close to 1.75, similar to that for both 
self-avoiding random walks and the backbone of a ran- 
dom percolation cluster in three dimensions . 

In Ref. J8|, it was shown that this quasi-one- 
dimensionality appears to arise from the tendency for 
mobile particles to follow one another. This is demon- 
strated in Fig. [H^, where we plot the time-dependent pair 
correlation function for the mobile particles, <?mm {r,t*) 
for different temperatures. At t = 0, this function co- 
incides with gjujvf(r) in Fig. O. For t > 0, the nearest 
neighbor peak moves toward r = 0, demonstrating that a 
mobile particle that at t = is a nearest neighbor of an- 
other mobile particle tends to move toward that particle 
at later times. We find that the peak at r = is highest 
near t = t*, and decreases for later times. A small but 
discernable peak at r — is also present in g(r,t*) Jtq] . 

Fig. [l7| shows a cluster of mobile particles at two dif- 
ferent times, t = and t — t* , to demonstrate the coop- 
erative, string-like nature of the particle motion. 

In a manner identical to our analysis of the mobile par- 
ticles, we define as immobile the 5% of the A particles 
that have the lowest value of [i. The pair correlation func- 
tion gn(r) between immobile particles shown in Fig. ^ 
shows that these particles also tend to be spatially corre- 



6 



lated. It is interesting to note that while the maxima in 
gii(r) are higher at all T than the corresponding max- 
ima in gAA the depth of the minima does not change 
appreciably for the lowest temperatures. Fig. [l9| shows 
the ratio T(r) — [gn(r)/ gAA^r)) — 1 as a function of r. 
In contrast to what we find for the most mobile parti- 
cles, the correlation between immobile particles does not 
show any evidence of singular behavior as T decreases. 
Instead, the correlation appears to grow and then "satu- 
rate" to some limiting behavior for all T < 0.468. More- 
over, Fig. [l9] shows that the local structure of the liquid 
appears to be more ordered in the vicinity of an immobile 
A particle than in the vicinity of a mobile A particle. 

In Fig. [2(] we show the size distribution of the clusters 
of immobile particles, formed with the same rule used 
for the mobile ones. One of the largest clusters found at 
T = 0.4510 is shown in Fig. ||. In the inset of Fig. |(j we 
show the mean cluster size S versus T — T c . We find that 
the mean cluster size of immobile particles is relatively 
constant with T. This may be because immobile parti- 
cles are relatively well-packed, and cannot grow beyond 
some limiting size ]47| . Or, these clusters may be the 
"cores" of larger clusters of particles with small displace- 
ments, that may grow with decreasing T, To elucidate 
this, more particles (e.g. the next 5% higher mobility) 
should be included in the analysis. We will return to this 
important point and provide further relevant data in the 
next section. 

The correlation between mobile and immobile parti- 
cles, measured by the pair correlation function gMi(f) 
(Fig. p2| ), shows that mobile and immobile particles 
are anti-correlated. A comparison between <?jvfj(r) and 
9aa(x)i shown in Fig. [2^, demonstrates that, over several 
interparticle distances, the probability to find an immo- 
bile A particle in the vicinity of a mobile one is lower 
than the probability to find a generic A particle. The 
figure also shows that the characteristic length scale of 
the anticor relation grows with decreasing T. This length 
scale does not show a tendency to diverge as T c is ap- 
proached. In particular, the curves for the two coldest 
runs (and closest to T c ) are almost coincident. 

VII. LOCAL ENERGY AND LOCAL 
COMPOSITION VS. MOBILITY 

We have seen in the previous section that despite the 
lack of a growing static correlation, a growing dynamical 
correlation — characterizing spatial correlations between 
particles of similar mobility — does exist. These corre- 
lations must therefore arise from subtle changes in the 
local environment that are not completely captured by 
the usual static pair correlation function. In this Sec- 
tion, we calculate several quantities to elucidate whether 
the mobility of a particle is related to its potential energy, 
and to the composition of its local neighborhood. 



In Fig. [24| we show the distributions of the potential 
energies of the 5% most mobile, 5% least mobile, and all 
particles at T — 0.4510, calculated at the beginning of an 
arbitrary time interval [t, t + t*]. The distributions have 
been normalized such that the area under each curve is 
one. The distributions differ by a small relative shift of 
the mean value, approximately 3% for the high mobil- 
ity distributions and somewhat less for the low mobility 
distribution. We find that the magnitude of the shift in- 
creases with decreasing T, but the relative shift appears 
to be independent of T . Since the liquid is in equilibrium, 
this shift will vanish for t — > 00. Thus, not suprisingly, 
mobile particles are those that in a time t* are able to 
rearrange their position so as to lower their potential en- 
ergy. It is worth noting that the mobility does not show 
any correlation with the kinetic energy of the particles 
measured at t = 0. The kinetic energy distributions of 
the subsets with different mobility coincide exactly with 
the average distribution, showing that the mobility can- 
not be related to the presence of "hot spots" in the liquid. 

We next divide the entire population of A particles into 
20 subsets, each composed of 5% of the particles. In the 
first subset we put the 5% of the particles with the highest 
values of fi (the mobile particles defined above), in the 
second subset the next 5%, and so on. The last subset 
thus contains the 5% most immobile particles. In Fig. |25| 
we plot (on the x-axis) the average mobility of each subset 
versus (on the y-axis) the average potential energy of that 
subset at t = 0. We find that the subset with the lowest 
mobility is also the one with the lowest potential energy. 
We also find that as the potential energy increases, the 
mobility increases. We see from the figure that the mobile 
particles are the subset with the highest average potential 
energy at t = 0. 

Two more points are worth noting in Fig. First, 
at all T the mobile particles move, on average, approx- 
imately one interparticle distance in the time interval 
[0, t*]. Second, for all T the difference in both mobil- 
ity and potential energy between the 5% most mobile 
particles and the next subset is significantly larger than 
between any other two consecutive subsets. This obser- 
vation suggests that the choice of 5%, while arbitrary, is 
a reasonable one. As shown in the figure, the separation 
between the 5% most mobile particles and the next sub- 
set shows a tendency to grow with decreasing T. Note 
however, that the distance between the lowest mobility 
subset and the next subset decreases with decreasing T, 
making it very difficult in the current approach to define 
an appropriate subset containing particles whose mobil- 
ity is distinctly lower than the rest. This, together with 
the result that the mean cluster size of immobile parti- 
cles is relatively constant over the range of temperatures 
studied, suggests that our analysis of the lowest subset 
is inadequate to fully characterize clusters of particles 
which do not move a substantial distance E§] . 

Thus we see that the gross structural information con- 
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tained in the potential energy is sufficient to establish a 
general correlation between energy and mobility. How- 
ever, as seen in Fig. |24], the distribution of potential en- 
ergies of mobile particles overlaps for most of the range 
of the abscissa with the distribution for the generic A 
particles. Thus, it is not possible to decide if a certain 
particle is mobile on the basis of that particle's instan- 
taneous potential energy alone. Other factors, such as 
defects in the local packing, and the relative potential 
energy of neighboring particles, must contribute as well. 

The relation between mobility and potential energy 
suggests a relation between mobility and local compo- 
sition. Indeed, a calculation of the pair correlation func- 
tions gMA{r) and <7mb(0 (see Figs. 26 - 12^) between a 



mobile particle and a generic A or B particle, respec- 
tively, shows that, on average, a mobile A particle tends 
to have less B particles, and more A particles, in its near- 
est neighbor shell than a generic A particle. 

A correlation can also be found between immobility 
and small composition fluctuations of the mixture. A 
comparison of <m_b(?") to the pair correlation function 
9ib{t), which measures the number of B particles a dis- 
tance r from a test immobile A particle, shows that an 
immobile A particle has, on average, more B particles in 
its nearest neighbor shell than does a generic A particle. 
As shown in Fig. |ol where the ratio [giB / gAB{r)] — 1 
is plotted as a function of r, this enhanced correlation is 
independent of T, and therefore does not suggest any ev- 
idence of A — B phase separation (recall that the chosen 
energy parameters preclude A — B phase separation). 

From these results it is clear how a correlation between 
mobility and local composition causes a correlation be- 
tween mobility and potential energy. Since the attractive 
interaction between A and B particles is stronger than 
either the attractive AA or BB interaction, the presence 
of a B between two A's reduces their potential energy. 
A particles in a B-rich region can thus be expected to 
have a reduced mobility. A particles in a B-poor region, 
however, will have a higher potential energy, resulting in 
a higher mobility. 



VIII. STRUCTURAL RELAXATION OF 
PARTICLE SUBSETS AND DYNAMICAL 
HETEROGENEITY 

We have shown that it is possible to select subsets of 
particles according to their maximum displacement over 
a timescale in the region of the late /3-early a relaxation. 
We have also shown that the particles belonging to sub- 
sets selected at the extrema of the mobility spectrum 
are spatially correlated, and are related to small fluctu- 
ations in the local potential energy, and, consequently, 
in the local composition of the mixture. All of the data 
presented here suggest that this supercooled liquid con- 
tains fluctuations in local mobility, with diffuse, quasi- 



one-dimensional regions of high mobility, and relatively 
compact regions of low mobility. 

To measure how long a mobile particle will continue 
to be mobile, we define a variable ^(t) as 1 if the i- 
th particle belongs to the subset of the 5% most mobile 
particles in the interval [t, t + t*], and otherwise. The 
function o~M{t) 



a M (t) 



1 



E 



Na 



(^(*k M (o)>-^), 0) 



measures the fraction of particles that are mobile in the 
interval [0,t*] and still mobile in the interval [t, t + t*], 
when the time origin is shifted by t. Here um is the 
number of mobile particles (320 in the present case), and 
Na is the total number of A particles (6400). The second 
term on the right-hand side of Eq. |9| is the number of 
particles that by random statistics would be classified 
as mobile in both time intervals. The normalization of 
<TM(t) is chosen so that ctm(0) = 1. The results for ctm(£) 
for the coldest T are shown in Fig. 52. 

We have also measured the fraction of particles that are 
immobile in the interval [0, t*] and still immobile in the 
interval [t, t + t*]. Analogous to the case for the mobile 
particles, we define o~i(t) as 



<7/(t) = 



1 



£<^H(o)> 



N A 



N A J> 



(10) 



where ni is the number of immobile particles (320) and 
v\ (t) is a function that is 1 if the i-th particle is an im- 
mobile one in the interval [t, t + t*], and otherwise. The 
function 07 (i) is also shown in Fig. |3l]. 

The functions o~M(t) and o~i(t) are memory functions of 
mobility. When they have decayed to zero, there are no 
particles that have retained memory of their mobility in 
the initial time interval. Because a particle's mobility is 
based upon a criterion that depends on t* , certain time- 
dependent functions measured for these subsets will have 
some "kink" at t* . If a different t* is chosen, the kink will 
move to the new t*. In this respect, there is no "natural" 
lifetime for these clusters — by definition, they survive 
for a time t* 

Nevertheless, we can obtain information from the form 
of the decay both before and after t* . Because the data 
was stored not less than every 3 time units, we are un- 
able to calculate ctm(£) and 07 (i) for t < 3. However, we 
see that these functions decay substantially before this 
time, since already at t = 3 both functions are signif- 
icantly smaller than one. After this initial short-time 
relaxation, a second decay of both functions is observed 
up to t = t* . At this time, a third decay process appears 
for the mobile particles, and possibly also for the immo- 
bile particles. The main point of Fig. |3l] is that beyond 
t* , both functions are less than 0.1. Thus there is only 
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a small tendency for particles to retain memory of their 
mobility beyond the initial time interval [5l| ] . 

Thus, mobile and immobile regions do not persist be- 
yond the time t* over which the particle mobility is mon- 
itored. After the observation time, mobile and immobile 
particles maintain little memory of their previous state. 
Therefore the strong correlations found between particles 
must arise from the motion itself. If, for instance, the mo- 
bility of a particle and its spatial correlation with other 
particles of similar mobility could be explained solely by 
local fluctuations in quantities like density or composi- 
tion, the mobility should persist until these fluctuations 
decay to zero. Instead, the dependence of the lifetime on 
the observation time can be explained if one assumes that 
particles can move only in a cooperative manner. Indeed, 
as was shown in Ref . M , clusters of mobile particles like 
that shown in Fig. |l5| can be decomposed into numerous, 
smaller string-like clusters ( "strings" ) of particles which 
follow one another in a cooperative fashion. 

Fig. |3^ shows the intermediate scattering functions 
F S M an d Fg(q,t), defined as the spatial Fourier 
transform of the self part of the van Hove correlation 
function G^ (r, t) or G], (r, t) of the mobile and immobile 
particles, respectively. Both functions are identical to 
the bulk F s (q,t) for times less than the "collision time" 
t = 0.09 The figure shows that a two-step relax- 

ation process occurs for the mobile particles, although 
the height of the plateau is smaller than for the bulk. 
The presence of the plateau in F S M (q, t) indicates that 
the mobile particles are subject to the same "cage effect" 
experienced by the other particles, although the effective 
cage "size" and "lifetime" are different. Thus clusters of 
mobile particles should not be thought of as "fluidized" 
regions of the liquid in the simple sense that those re- 
gions might behave like high temperature or low density 
liquids. Instead, the difference between the mobile par- 
ticles and the rest of the sample appears to be, from the 
point of view of the single particle dynamics, that they 
"escape" the cage earlier than the other particles. 

We also see that the three curves in Fig. ^ cannot be 
superimposed by scaling the time axis in the same way as 
one can superimpose F(q,t) curves for different temper- 
atures. Again, this indicates that the mobile and immo- 
bile subsets are not simply "hotter" or "colder" subsets of 
the sample, in agreement with the perfect superposition 
of the kinetic energy distributions. 

In contrast to the bulk average F s (q, t), F^ 4 (q, t) is not 
a monotonically decreasing function of time. For times 
longer than t*, a small but clearly detectable increase of 
the function can be noticed in Fig. ^2|. This behavior can 
be interpreted as a tendency of a small fraction of the 
particles that we have selected to return towards their 
position at the beginning of the selection interval. These 
particles may also be those that contribute to the small 
memory effect observed in <JM{t) in Fig. but further 
analysis is required to establish this connection. 



Finally, we show in Fig. |33| the fraction <fi of particles 
that at time t have not yet been labeled mobile. This 
function is calculated by labeling the mobile particles in 
the first interval [0, f ], and then shifting the interval by 
t and reassigning the particle mobilities. Thus at t = 0, 
95% of the particles have not been labeled mobile. In the 
interval [t,t + t*], more particles will have been labeled 
mobile, so 4> will decrease. We have normalized <f>(i) such 
that 0(0) = 1. Also shown in Fig. [53| is the long-time 
a-relaxation part of the bulk F s (q,t) for q = q max - Fits 
to both functions are also shown. Both functions fit well 
to a stretched exponential y(t) = Aexp[(— t/r)P] with 
(3 = 0.75 and r Q = 655 for the intermediate scattering 
function, and with (3 = 0.78 and r a = 475 for <fi. That 
both functions have a similar form (similar /3), and sim- 
ilar time constants, suggests that the process by which 
immobile particles become mobile governs the long time 
structural relaxation of density fluctuations at wavevec- 
tors corresponding to the peak of the static structure 
factor. Moreover, it demonstrates that the "arbitrary" 
choice of 5% represents a physically meaningful fraction 
of the system . 

IX. DISCUSSION 

In this paper, we have described an investigation of the 
individual particle dynamics of a cold, dense Lennard- 
Jones mixture well above the glass transition in an effort 
to discover if the liquid is dynamically heterogenous, and 
if so to determine the extent and nature of the dynami- 
cal heterogeneity. Since there were no quantitative the- 
oretical predictions regarding this matter, the approach 
we have taken is exploratory; particle trajectories were 
saved during the course of the simulation and then an- 
alyzed and visualized in numerous ways. We find that 
this supercooled liquid is "dynamically heterogeneous" 
because particles with similar mobility are spatially cor- 
related. Note that our definition of heterogeneity is dif- 
ferent from the one used, for instance, in 4-D NMR ex- 
periments, where the system is defined as heterogeneous 
if a slow subset remains slow for times longer than the 
average relaxation time Qj . We further find that highly 
ramified clusters of mobile particles grow with decreasing 
T and appear to percolate at the mode-coupling temper- 
ature. This is the first evidence for a percolation transi- 
tion coincident with T c , and it is very different from the 
type of percolation transition proposed in free volume 
theory [pa]. It is especially interesting since MCT docs 
not make any predictions regarding clusters or diverging 
length scales. We also find that particles of low mobil- 
ity form relatively well-ordered, compact clusters which 
do not appear to grow with decreasing T if the num- 
ber of immobile particles included in the subset is kept 
constant. Although mobile and immobile clusters are 
anti-correlated, there is no tendency towards bulk phase 
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separation of mobile and immobile regions because of the 
highly ramified, extended nature of the mobile regions. 

In our analysis, we find no evidence to support a pic- 
ture in which the system can be thought of as a collection 
of subvolumes that each relax independently and simul- 
taneously with their own time constant. Instead, it ap- 
pears that at any given time, most particles are localized 
in cages and a small percentage of particles form large 
clusters of smaller, cooperatively rearranging "strings." 
After rearranging, these mobile particles become caged 
themselves, and others become mobile. This process re- 
peats until, on the time scale of the a relaxation, each 
particle has rearranged at least once. Thus the structural 
relaxation of the liquid appears to be highly cooperative 
in the spirit of Adam and Gibbs, but where different sub- 
volumes of the liquid are able to relax only after other 
subvolumes relax. This will be further explored in a sep- 
arate publication Q . 
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FIG. 1. Incoherent (self) intermediate scattering function 
F s (q,t) for q max = 7.12. 
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FIG. 2. Pair correlation function gAA(r) of the A particles 
for three different temperatures. 
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FIG. 3. Pair correlation function gAB(r) between A and B 
particles for three different temperatures. 



FIG. 5. Mean square displacement (r (t)) of the A particles 
vs. time. 
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FIG. 4. Pair correlation function pss(r) of the B particles f 
for three different temperatures. 

FIG. 6. Solid line: 47rr 2 G s (r, t) of the A particles for three 
times at T = 0.4510. Dashed line: Gaussian approxima- 
tion calculated using the measured (r 2 (t)) for the same three 
times. 
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FIG. 7. Same intermediate time data as in previous figure, 
but enlarged. Solid line: 4irr 2 G s (r, t) of the A particles at 
t — 155.5 at T = 0.4510. Dashed line: Gaussian approx- 
imation calculated using the measured (r 2 (t)} for the same 
time. 




FIG. 9. Probability distribution An/^Pip, t*) (dashed line) 
of a particle having a maximum displacement of magnitude fj, 
at t*. For comparison, the distribution 4irr 2 G s (r, t*) is also 
shown. 
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FIG. 11. Pair correlation function gMM(r) between mobile 
A particles at four different temperatures. 
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FIG. 13. Distribution of the size n of clusters of mobile 
particles. Inset: Mean cluster size S plotted versus T — T c , 
where is the fitted MCT critical temperature T c — 0.435. The 
straight line is a power law fit S ~ (T — T c )~"' , with 7 = 0.618. 
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FIG. 12. T(r) = [gMMir)/g A Air)} - 1 vs. r for different 
temperatures. 
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FIG. 14. Mean cluster size S plotted versus T, for subsets 
containing 3%, 5% and 7% of the most mobile particles. The 
data for the 5% are the same as those shown in the inset of the 
previous figure. The lines are power law fits S ~ (T — T v ) ' . 
Best fit parameters are T p — 0.440, 0.431 and 0.428, respec- 
tively, and 7 = 0.397, 0.687, and 0.741, respectively. 
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FIG. 15. One of the largest clusters of mobile A particles 
found at T=0.4510. The cluster is composed of 125 particles, 
which are represented here as spheres of radius r = 0.5<7 a a. 
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FIG. 18. Pair correlation function gn(r) between immobile 
particles. 
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FIG. 19. T(r) = [gn(r)/g A A(r)] - 1 vs. r for different 
temperatures. 



FIG. 16. Time-dependent pair correlation function 
gMM{r,t*) vs. r, for different temperatures. 



FIG. 17. A cluster of mobile particles at t = (light 
spheres) and t = t* (dark spheres), for T = 0.4510. 
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FIG. 20. Distribution of the size n of clusters of immobile FIG. 22. Pair correlation function g M i{r) between mobile 

particles. Inset: mean cluster size S plotted vs. T - T c . an d immobile A particles. 
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FIG. 24. Distribution of the potential energy of all the A 
particles, of the mobile A particles and of the immobile A 
particles for T = 0.4510. 



FIG. 26. Pair correlation function gMA(r) between mobile 
A and generic A particles. 
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FIG. 25. Potential energy (Z7j) as a function of the mobility 
for the A particles. The A particles have been divided 
into 20 subsets according to their mobility at t* . Each subset 
is represented by a point in the graph. The energy scale for 
T — 0.4510 is on the left hand side y axis, while the energy 
scale for T = 0.550 is on the right hand side y axis. 
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FIG. 27. Pair correlation function gMs(r) between mobile 
A and generic B particles. 



17 




-0.1 



T = 0.5495 
T = 0.5052 
T = 0.4795 
T = 0.4572 
T = 0.4510 



6 



FIG. 28. T(r) = [g M A/gAA(r)] - 1 vs. r for different tem- 
peratures. 
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FIG. 30. T(r) = [gi B /gAB(r)] - 1 vs. r for different tem- 
peratures. 
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FIG. 29. T(r) = [g M B/gAB(r)} - 1 vs. r for different tem- 
peratures. 
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FIG. 31. <JM{t) (dot-dashed line) and ai(i) (solid line) for 
T=0.4510. 
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FIG. 32. Self intermediate scattering function for 
T=0.4510 for all the A particles (solid line), for the mo- 
bile particles (dashed line) and for the immobile particles 
(dot-dashed line). 
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FIG. 33. 4>{t) (dashed line) and F s (q,t) (dotted line) for 
T=0.4510. The solid lines are fits to a stretched exponential. 
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